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Abstract. Consider the (simplified) Leslie-Erickson model for the flow of nematic liquid crystals 
in a bounded domain Q C M n for n > 1. This article develops a complete dynamic theory for these 
equations, analyzing the system as a quasilinear parabolic evolution equation in an L p — L^-setting. 
First, the existence of a unique local strong solution is proved. This solution extends to a global 
strong solution, provided the initial data are close to an equilibrium or the solution is eventually 
bounded in the natural norm of the underlying state space. In this case the solution converges 
exponentially to an equilibrium. Moreover, the solution is shown to be real analytic, jointly in time 
and space. 



1. Introduction 

We consider the following system modeling the flow of nematic liquid crystal materials in a bounded 
domain SI C W l 

d t u + (it • V)u - j/Am + Vtt = -Adiv([Vd] T Vd) in(0,T)xO, 

d t d+{u-V)d = j(Ad+\Vd\ 2 d) in(0,T)xfi, 

(1.1) { divu = in(0,T)xfi, 

(u,d v d) = (0,0) on (0,T)x9O, 

(u,d)\ t =a = (u ,d ) in SI. 

Here, the function u : (0, oo) x £1 — > W 1 describes the velocity field, 7r : (0, oo) x £1 — » K the pressure, 
and d : (0, oo) X CI — >■ 1" represents the macroscopic molecular orientation of the liquid crystal. Due 
to the physical interpretation of d it is natural to impose the condition 

(1.2) \d\ = l in(0,T)xn. 

We will show in the following that this condition is indeed preserved by the above system; see 
Proposition 14.31 below for details. 

The constants v > 0, A > 0, and 7 > represent viscosity, the competition between kinetic energy 
and potential energy and the microscopic elastic relaxation time for the molecular orientation field, 
respectively. For simplicity, we set v = A = 7 = 1, which does not change our analysis. 

The continuum theory of liquid crystals was developed by Ericksen and Leslie during the 1950's 
and 1960's in j9j [17] • The Ericksen-Leslie theory is widely used as a model for the flow of liquid 
crystals, see for example the survey articles by Leslie in [TU] and also [H O Q31 [20] • 

The set of equations (|1.1[) was considered first in [21] , however for the situation where in the second 
equation of f| 1 . 1 [) the term |V<i| 2 ci is replaced by f(d) = \7F(d), i.e. 

d t + (u-V)d = j(Ad-f(d)), 
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where F : K 3 — > K is a smooth, bounded function. Note that in this situation, the condition (|1.2[) 
cannot be preserved in general. Thus, this condition was replaced in |20| and |21| by the Ginzburg- 
Landau energy functional, i.e. / is assumed to satisfy f(d) = VF(d) = V^(\d\ 2 - l) 2 . In 1995, 
Lin and Liu |21) proved the existence of global weak solutions to (jl.ip in dimension 2 or 3 under the 
assumptions that uq G L2(fi),do G and do G H 3 ^ 2 (dfl). Existence and uniqueness of global 

classical solutions was also obtained by them in dimension 2 provided uq G do G H 2 (fl), and 

provided the viscosity v is large in dimension 3. For regularity results of weak solutions in the spirit 
of Caffarelli-Kohn-Nirenberg we refer to [22] . 

Hu and Wang [H] considered in 2010 the case of f(d) = and proved existence and uniqueness 
of a global strong solution for small initial data in this case. They proved moreover that whenever 
a strong solutions exist, all global weak solutions as constructed in [21 must be equal to this strong 
solution. The idea of their approach was to consider the above system (jl.l[) as a semilinear equation 
with a forcing term Adiv([V<i] T V<i) on the right-hand side. 

The full system (|1.1[) with f(d) = \\7d\ 2 d was revisited by Lin, Lin, and Wang in 2010. They 
proved in [19] interior and boundary regularity theorems under smallness condition in dimension 2 
and established the existence of global weak solutions on bounded smooth domains O C I 2 that 
are smooth away from a finite set. Furthermore, Wang proved in [21] global well-posedness for this 
system for initial data being small in BMO^ 1 x BMO in the case of a whole space, i.e. ft — R", by 
combining techniques of Koch and Tataru with methods from harmonic maps to certain Riemannian 
manifolds. 

For results on the compressible case we refer to [TSl [53]. Here, local existence of strong solutions 
is proved. The latter turn out to be even local classical solution. 

Summarizing, we observe that in particular results for local as well as global, strong solutions in 
the three dimensional setting for the full system (jl.ip . obeying also the condition (jl.2[) . do not seem 
to exist so far. 

Recently, Li and Wang claimed in [18] such a result. More precisely, they claimed the existence 
and uniqueness of a strong solutions to (|1.1[) in bounded, smooth domains (however, not satisfying 
(jl.2p ). Their idea was to rewrite (jl.ip as a semilinear equation for the Stokes equation coupled to 
the heat equation with a right hand side of the form 

F{u, d) := (-(u • V)it - div([Vd] T Vd), -{u ■ V)d + \Vd\ 2 d) . 

Unfortunately, their approach and their main result [181 Theorem2.1] relies on an incorrect regularity 
property for the solution of the heat equation [TBI Theorem 3.1]. This result would imply further 
regularity properties for d and hence for F(u, d), which however are not true. Note that the (incorrect) 
assertion of [TSl Theorem 3.1] is crucial for their approach. Thus, the theory for local as well as for 
global strong solutions to (jl.ip . also satisfying (|1.2p . needs clarification. 

It is the aim of this paper to present a complete theory for global strong solutions to (jl.ip satisfying 
(|1.2D as well as for their dynamical behaviour in the n-dimensional setting, where n > 1. 

Our main idea is to consider (jl.ip not as a semilinear equation as done in all of the previous 
approaches but as a quasilinear evolution equation. We thus incorporate the term div([Vd] T V<i) into 
the quasilinear operator A given by 



A(d) = 



A PB{d) 
V 



where A denotes the Stokes operator, V the Neumann Laplacian, and B is given by 

[B{d)h]i := didiAhi + dkdidkdihi, 

for which we employ the sum convention. Note that B(d)d = div ([\7d] T \7d). 

We then develop a complete dynamic theory for (jl.ip - (|1.2p . In fact, first by local existence theory 
for abstract quasilinear parabolic problems, we prove the existence and uniqueness of a strong solution 
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to l|l.l|) - (|1.2j) on a maximal time interval. Thus, (|l.ip - (|1.2p give rise to a local semi-flow in the natural 
state space. 

Furthermore, the equilibria E of (|l.ip - (|1.2p are determined to be 

£ = {(0,4) : 4 el", \d*\ = 1}, 
and the well-known energy functional 

E = i / [\u\ 2 + \Vdf}dx 

for (|l.ip - (|1.2p is shown to be a strict Ljapunov- functional. In addition, the equilibria are shown to be 
normally stable, i.e. for an initial value close to £, the solution of (|l.ip ~ (|1.2p exists globally and the 
solution converges exponentially to an equilibrium. More generally, a solution, eventually bounded 
on its maximal interval of existence, exists globally and converges to an equilibrium exponentially 
fast. 

Due to the polynomial character of the nonlinearities, we can even show that the solution of 
(|l.ip - (|1.2p is real analytic, jointly in time and space. 

Our approach is based on the theory of quasilinear parabolic problems and relies in particular on 
the maximal L p -regularity property for the heat and the Stokes equation. In particular, we refer here 

to [U 00 mug EH m EH]. 

The plan for this paper is as follows. We begin by collecting general results from the theory of 
quasilinear parabolic evolution equations. Then, in Section [3] we introduce our formulation of (jl.ip . 
Section @] deals with local well-posedness and regularity of solutions to (|l.ip - (|1.2|) ; in particular we see 
that the solution is real analytic. The generalized principle of linearized stability yields the stability of 
equilibria and convergence of solutions is proved in Section [5l Moreover, by means of the associated 
energy functional, we prove convergence of a solution to an equilibrium, whenever the solution is 
eventually bounded in the natural state space. 

2. Quasilinear Evolution Equations 

Let Xq and X\ be Banach spaces such that X\ Xq, i.e. X± is continuously and densely embedded 
in Xq. Let J — [0, a] for an a > 0. By a quasilinear autonomous parabolic evolution equation we 
understand an equation of the form 

(QL) z(t) + A(z(t))z(t) =F(z(t)), t€J, z{0) = z o , 

where A is a mapping from a real interpolation space X 14l with suitable weights between Xq and X\ 
into C{Xq, X\). Our approach relies on the maximal L p -regularity of A(v) for v G X y ^. For details 
we refer e.g. to [?]• 

The equation ( |QL[ ) is investigated in spaces of the form L p (J;Xq) with temporal weights. More 
precisely, for p G (1, oo) and /j, G (1/p, 1], the spaces L p ^ and H p are defined by 

L v ^{J;Xx) := {z: J -> X x : t^z G L^J-X^}, 

HIJJ-Xq) := {z G L p J J; X )n Wl(J;X ): z G L PiP ( J; X )}. 

It is clear, that 

L P (J;X) L Ptli (J;X) and L p ([0, a]; X) ^ L Ptfl ([r, a]; X), 

for all Banach spaces X and r G (0, a). It has been shown in [26l Theorem 2.4] that L p -maximal 
regularity implies also Lp !M -maximal regularity, provided p G (1, oo) and fi G (1/p, 1]. The trace space 
of the maximal regularity class containing temporal weights, 

z&H^(J;X )nL p>M (J;Xi) 
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has been characterized in |26[ Theorem 2.4] as 

Xy,H = (Xq, -^l)ft-l/p,p) 

provided p e (1, oo) and p G (1/p, 1]; see also [2H Theorem 4.2]. 
We now impose precise assumptions on A and F. 

(A) A £ C w (X r /1 ; £(JTq, Xi)), and has maximal L p -regularity for each v G X 7iM . 
(F) i^eC w (X T ' M ;Xo). 

Even under less restrictive Lipschitz type assumptions on A and F, local in time existence of ( |QL[ ) 
was shown by Clement-Li [5] in the case = 1 and by Kohne-Priiss-Wilke [16, Theorem 2.1, Corollary 
2.2] for the case G (l/f>, 1]. 

Proposition 2.1. Le£ 1 < p < oo, /i G (1/p, 1], z G X^ t a, and suppose that the assumptions (A) and 
(F) are satisfied. Then, there exists a > 0, such that ( |QL[ ) admits a unique solution z on J = [0, a] 
m the regularity class 

z G Hi „( J; X ) n L P , M (J; Xi) C(J; Jf 7iM ) n C((0, a]; X 7 ). 

T/ie solution depends continuously on zq, and can be extended to a maximal interval of existence 
J(z ) = [0,t+(z )). 

Parabolic problems allow for additional smoothing effects. In this respect, a method due to An- 
genent [3] is well known. We only state here a variant of it which is adapted to ( |QL[ ) ; see [25l Theorem 
5.1] for the case \i = 1. By a slight adjustment of its proof to the situation of temporal weights, this 
result remains true also for maximal regularity classes using this type of weights. 

Proposition 2.2. Let 1 < p < oo, fi G (l/p,l], a > 0, and assume that (A) and (F) hold. Let 
z G Hp Xq) H L p . M (J; Xi) be a solution of ( |QL[ ) on J — [0,a]. Then 

t k [±] k z G HlJJ;X ) n L^JiXt), k G N. 
Furthermore, z is real analytic with values in X\ on (0, a). 
We denote the set of equilibria of ( |QL[ ) by 

£ = {z* G Xi : A(z*)z r — F(z*)}. 

The following result on global existence and stability of was proved in [281 Theorem 2.1] assuming 
only that A and F are of class C 1 . 

Proposition 2.3. Let 1 < p < oo and assume that assumptions (A) and (F) hold. Let A$ be the 
linearization of ( |QL[ ), i.e. let 

Aqw = A(z*)w + (A (u*)w)u* — F'(u*)w, w e Xi. 

Suppose that G £ is normally stable equilibrium, i.e. 

(i) near u* the set of equilibria £ C X% is a C 1 -manifold in X± of dimension m G No, 

(ii) the tangent space of £ at is given by N(Aq), 

(hi) is semi-simple eigenvalue of Aq, i.e. N(Aq) © R(Aq) — Xq, 
(iv) a(A ) \ {0} C C+ = {x G C: Rca; > 0}. 

Then w* is stable in X^ . Further, there exists a number p > such that the unique solution z of ( |QL[ ) 
with initial value z G By (0, p) exists on R + and converges at an exponential rate to some m^. G £ 
in X 7 as t oo. 



We finish the section with another result on global existence result for (QL|; see [16, Theorem 3.1]. 
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Proposition 2.4. Let 1 < p < oo, ji 6 (1/p, 1], zo G -X-y,|» arl ^ ^ ^ = [0>o] or J = K+. Suppose 
that assumptions (A) and (F) are satisfied and that the embedding X 1 ^ <—> X 1 is compact. Assume 
furthermore that the solution z of ( |QL[ ) is eventually bounded in X 1 on its maximal interval of 
existence, i.e. that z satisfies 

zeBCIM+tzo));!,) 

for some r 6 (0,t + (zo))- Then the solution z exists globally and for each 8 > 0, the orbit {z(t)} t >s 
is relatively compact in X 7 . If in addition zq G Xj, then {z(t)} t >o is relatively compact in X 1 . 

3. Nematic liquid crystals AS QUASILINEAR EVOLUTION EQUATIONS 

We now reformulate (jl.ip equivalently as a quasilinear parabolic evolution equation for the un- 
known z = (u, d). To this end, for 1 < q < 00 define the Banach spaces X by 

X := L, i<r (fl) x L q (Q) n , 

where C P n is a bounded domain with boundary dil S C 2 . The subscript er in L qrT (Q) means as 
usual the subspace of L q {Q) n consisting of solenoidal vector fields. 

The Neumann-Laplacian T> q in L q (Q) is defined by T> q = —A with domain 

D(V q ) := {d e H 2 q {ti) n : d„d = on SQ}. 

It is well-known that V q has the property of L p -maximal regularity; see pj Theorem 8.2]. 

Let P : L q (Ct) n — > L g)CT (f2) denote the Helmholtz projection. We then define the Stokes Operator 
A q = —PA in L QiCr (Cl) with domain 

D(Aq) = {u e H 2 q {VL) n : divu = in Q, u = on dtt}. 

It is also well-known that A q has the property of L p -maximal regularity; see e.g. [29) , IT2 | [TT]. 
Next, we define the space Xi by 

Xx := D{A q ) x D(V q ), 
equipped with its canonical norms. Then X\ Xq densely. 



The quasilinear part A(z) of ( |QL[ ) is given by the tri-diagonal matrix 

Aq PB q (d) 



A{z 

where the operator B q is given by 



[B q (d)h]i := didiAhi + d k did k dih h 
for which we employed the sum convention. Note that 

= div([Vd] T Vd). 

Obviously, S(d) : X\ — \ Xq is bounded for each d S C 1 (f7)™ and the map d M- PB g (d) is polynomial, 
hence real analytic. By the tri-diagonal structure of A(z) and by the regularity of B one can easily see 
that A{z) also has the property of L p -maximal regularity, for each z e C 1 (il) 2 ". Indeed, for a fixed 
right-hand side (f u ,fd) G ^p(0, a; X 7 M ) and initial values (ito,<io) G X^^, we may use the maximal 
regularity of T> q to obtain a solution d of the heat equation with Neumann boundary condition in the 
right maximal regularity class. By setting 

J u := f u -PB q (d)d 

as right-hand side for the Stokes equation, we obtain a solution u in the right maximal regularity 
class due to the fact that B q (d) is linear and bounded. 
The right-hand side F(z) of ( |QL[ ) is defined by 

F(z) = (-Pit • Vit, -u • Vd + I Vd| 2 d), 
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which is also polynomial, hence a real analytic mapping from C 1 (^l) 2n into Xq. 
Note that (A) and (F) hold, as soon as we have the embedding 

The space A 7 is given by 

X 1 := (X ,Zi)i_i/ PiP = D Aq {\ - l/p,p) x D Vq (l - l/p,p); 

see pQ[S]. As explained in Section [3J we consider L p -spaces with temporal weights. The trace space 
of the class 

zGi7 p 1 jM (J;A )nL p , M (J;X 1 ) 

now reads 

X 1 , f t = (X Q , Xi) M _i/p, p = D Aq (fi - l/p,p) x D Vq {n - l/p,p), 
provided p G (1, oo) and fi £ (1/p, 1]; see [24j Theorem 4.12]. 

In order to obtain the embeddings X 1 C 1 (Vt) 2n and more generally A 7iM <—t C 1 (il) 2n we impose 
on p, q £ (1, oo) now the conditions 

2 n 1 1 n 

(3.1) - + -<1, - + - + — < M <1. 

P 9 2 p 2q 

Standard Sobolev embedding theorems can then be applied. 

Further, we recall from [301 Theorem 4.3.3] and [21 Theorem 3.4], respectively, the following char- 
acterizations of the interpolation spaces involved, 

d G Dd,(a* - l/p,p) ^ d £ B 2 £- 2/p (n) n , d v d = on 90, 

and 

u G D Aq Qi - l/p,p) & u £ B^- 2 /f(n)« n L 9 , CT (0), u = on 90. 

Observe that both of these characterizations make sense, since the condition (|3.1[) guarantees the 
existence of the trace. 

4. Existence, uniqueness, and regularity of solutions 

We start this section by applying Proposition 12.11 to obtain the following result on local well- 
posedness of 

Theorem 4.1. Letp,q,fj, be subject to (|3.1| . and assume zq = (uo,e?o) £ -^t,mj w hich means that 
UQ,d £ B 2 p 2 ^ p (0)™ satisfy the compatibility conditions 

divwo = in O, wo = d u do = on 90. 

T/ien for some a = a(zo) > 0, there is a unique solution 

z£H^(J,X )nL p ^(J;X 1 ), J=[0,a], 

of (|1.1| on J. Moreover, 

z e C([0,a]]Xw)nC((p,a];X y ), 
i.e. the solution regularizes instantly in time. It depends continuously on zq and exists on a maximal 
time interval J(z ) = [0,t + (z o )). Therefore problem (|l.ip , i.e. ( |QL[ ) 7 generates a local semi-flow in 
its natural state space X^ ^. 

Remark 4.2. Assuming that 2/p + n/q < 1, for e > we may choose /x subject to (|3.1I) such that 

H 1 q + ^ +£ {n) n ^ Bf;- 2 / p (n) n H l q + ^~ e {ny i 

due to Sobolev embeddings [30 , Theorem 4.6.1]. Furthermore, we can choose p, q large with 

c l+£ {n) n ^ B 2 %- 2 /p(n) n . 
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Employing different time weights for u and d, an inspection of the above proofs shows that the 
assertion of the above theorem remains true provided uo G C a (Q). 

The following result tells that the condition (|1.2[) is preserved by 

Proposition 4.3. Suppose that n,p,q are satisfying (13. ip and let z = (m ,c?o) G -^y.u uwift |dol = 1> 
a > 0. Lei 

zG^(J;X )nL p , M (J;X 1 ) 

6e a solution of on i/ie interval J = [0, a]. TTien |d(t)| = 1 ftoZds /or a// i G [0, a]. 

Proof. Setting 93 = |d| 2 — 1 the elementary identities, 

d t \d\ 2 = 2d ■ d t d, A\d\ 2 = 2Ad • d + 2|Vd| 2 , V|d| 2 = 2d-Vd, 

and multiplication with d of the second line in Ijl.ljl yields the problem 

dtf + u-Vif = Aip + 2\Vd\ 2 ip mil 

d v ip = on 9f2, 

ip{0) = in ft, 

provided | rfo I = 1- Uniqueness of this parabolic convection- reaction diffusion equations yields <p = 0, 
i.e. \d\ = 1. □ 

As the nonlinearities A and F are real analytic we may employ Angenent's method (Proposition ^. 2p 
to obtain further regularity of the solutions of (|l.ip . 

Proposition 4.4. Suppose that (J>,p,q satisfy (|3.ip . zo S X-y,n> anc ^ a > and /ei 

ze^(J;A )nL p , M (J;A 1 ) 
be a solution of on i/ie interval J — [0, a]. TTien /or eac/i G N, 

i fc [^] fc z G Hl^{J;X ) n L p , M (J; X x ). 

Moreover, z G C w ((0, a); Xl). 

We will employ Proposition 14.41 in the following to justify the regularity of time derivatives of the 
energy functional. 

Remark 4.5. Employing scaling techniques jointly in time and space, it is possible to show via 
maximal regularity and the implicit function theorem that u, tt, d are real analytic in (0, t + (zo)) x 57; 
see |251 Section 5] for parabolic problems, and specifically for a Navier-Stokes problem [27]. As we 
will not use this result below we omit the details, here. 

5. Stability and Convergence to Equilibria 

We consider the set £q — {0} x R n , which are obviously equilibria of (II. ip . This set forms a 
n-dimensional subspace of X\, hence a C 1 -manifold with tangent space {0} x K™ at each point 
(0, o?*) G £0 ■ The linearization of (jl.ip at z* G £0 is given by the linear evolution equation 

z + A*z = /, z(0) = zq, 

in Xq, where 

A* = dmg{A q ,V q ), D(A*) = X v 

As Q is bounded, the spectrum o~{A q ) consists only of positive eigenvalues and (jL o~(A q ). On the 
other hand, T> q has as an eigenvalue, which is semi-simple and the remaining part of o~(T> q ) consist 
only of positive eigenvalues. Thus o~(A*) \ {0} C [S, 00) for some 5 > and the kernel of A* is given 

by 

N(A*) = {0} x M™, 
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which equals the tangent space. In Remark 2.2 in [28J it is shown that all equilibria close to 2* are 
contained in a manifold A4 of dimension n = dim(7V (A*)) . Since the dimension of £0 is also n, there 
exists an open set V C X\ with Mfl7 = ^ny = £onl''; i.e. 8 H V contains no other equilibrium. 
As a result we see that the equilibrium is normally stable. 

Now we are in position to apply Proposition 12.31 to conclude the following stability result for the 
equilibria of 

Theorem 5.1. Let p,q satisfy the first inequality in (|3.1ll . Then each equilibrium z„ G {0} x M™ is 
stable in Xy, i.e. there exists e > such that a solution z(t) of (jl.ll) with initial value zq G X^, 
\zq — Z*| y < e, exists globally and converges exponentially to some z^ G {0} x R" in 

We next consider the energy of the system given by 
(5.1) E = ~ ^\\u\ 2 + \Vd\ 2 ]dx = E kin + E pot . 

Let z — Hp ( J; Xq) n L P , M ( J; X\) be a solution of f| 1 . 1 [) (| 1 . 2(1 on J = [0, a]; according to Proposition 
4.4 it belongs to C 1 ((0, a); Xi). Using sum convention we have by an integration by parts 

^-E kln {t) = [ d t u ■ udx 
at J n 

= / [-(u ■ V)u - Vtt + Ait - div([Vd] T Vd]) ■ udx 
Jn 

= — I \\7u\ 2 dx+ / dkdjdidjdkUidx, 
Jn Jn 



as div u = in Q and u = 0on dfl. On the other hand, we have by another integration by parts 

f \Ad+\Vd\ 2 d\ 2 dx= [ [Ad+\Vd\ 2 d][d t d+(u-V)d]dx 
Jn Jn 

= - [ [dtVd : Vd- \\7d\ 2 d t \d\ 2 /2]dx 
Jn 

+ [ [(u- V)d- Ad+\\7d\ 2 {u- V)\d\ 2 /2}dx 
Jn 



Ep t(t) - / dk(uididj)dkdjdx 



dt 



= ~4;Epot(*) - / dkUtdidjdkdjdx, 
dt Jq 



by |d| = 1 and the Neumann boundary condition for d. Combining these equations, we obtain the 
energy identity 

d f 

(5.2) —E(t) = - [\Vu\ 2 + \Ad+\Vd\ 2 d\ 2 }dx. 

dt Jq 

Therefore E(i) is non-increasing along solutions. But E is also a strict Ljapunov functional, i.e. strictly 
decreasing along non-constant solutions. In fact, if dE(t)/dt — at some time instant, then by the 
energy equality we have Vu — and Ad+ \ Vd\ 2 d — in fi. Therefore u = by the no-slip condition 
on dQ, and d satisfies the nonlinear eigenvalue problem 

C Ad+\Vd\ 2 d = infi, 

(5.3) I \d\ 2 = 1 in f2, 

[ dud = on <9f2. 

But, as the lemma below shows, this implies Vd = in f2, hence d = d* is constant and := (0, d*), 
|d*| = 1 is an equilibrium of the problem. 
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Lemma 5.2. Suppose that d G iLf(f2;R") satisfies (|5.3p . Then d is constant in ft. 

Proof. The idea is to reduce inductively the dimension N = n of the vector d. This can be achieved 
by introducing polar coordinates according to 

d\ = ci cos 8, d2 — ci sin#, dj — Cj-i, j > 3. 

Simple computations yield 

l = M 2 Hc| 2 , |Vd| 2 = |Vc| 2 +c 2 |V0| 2 , 

and 

A Cj + [|Vc| 2 + c 2 |V6f] Cj = in ft, 
as well as d v Cj = on dfl for j = 2, . . . , n — 1. Moreover, by an easy calculations we further obtain 

-Acx + Cl |V6»| 2 = [|Vc| 2 + c 2 |V6»| 2 ] Cl in O, 

and 

ciA0 + 2Vci • = in n, 

as well as 

S^Ci = c\d v B = on 90. 
Multiplying the former equation by C\Q and integrating over fi we deduce 

= / [ciA0 + 2Vci • V6)]a6dx = / div [c 2 V6>]6>dx = - / c\\X/0\ 2 dx, 
Jn Jn Jn 

hence ciV# = 0. This implies that c satisfies the problem (|5.3[) where the vector c has dimension 
iV — 1. Inductively, we arrive at dimension N = 1 and if d is a solution of (|5.3p with dimension 1, 
then d = 1 or d = — lby connectedness of f2. □ 

Note that the side condition \d\ = 1 is important at this point. Summarizing we proved the 
following result. 

Proposition 5.3. The energy functional E defined on X 1 is a strict Ljapunov function for system 
(|l.ip - (ll.2p . The equilibria of this system are given by the set 

£ = {z* = («*, d*) : u* = 0, d* G M™, |d*| = 1}, 

which forms a manifold of dimension n — 1 . The corresponding pressures are constant as well. 

Suppose finally that z is a solution of (|1.1|) - (|1.2|) which is eventually bounded in X 1 on its maximal 
interval of existence. Then, by Proposition 12.31 this solution is global and z([S, oo)) C X 7 is relatively 
compact. Therefore its limit set 

lu(zo) = {«£ X~f : 3t n j" oo s.t. z(t n ; zo) — > v in X 7 } 

is nonempty. As E is a strict Ljapunov functional for (|1.1[) - (|1.2[1 . we obtain dist(z(t, zq),lu(zo)) —> 
in Xy for t —> oo and lu(zq) C £ C Xi. Now Theorem 15.11 applies and we may conclude that 
— > Zoo G £ in X 7 as t — > oo. In summary we proved the following result. 

Theorem 5.4. Lei fi,p,q satisfy (|3.ip . Lei zq = (ito,do) G wi/i |do| = 1 and suppose that the 

solution z(t) of (jl.ip is eventually bounded in X 1 on its maximal interval of existence, i.e. 

z G BC{[r,t+{z ))-X 1 ) 

for some r G (0,i + (zo)). Then z{t) exists globally and converges to an equilibrium z^ G £ in X 1 , as 
t — > oo. The converse is also true. 
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